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Resonance modes in single crystal sapphire (a-A^Os) exhibit extremely high electrical and me- 
chanical Q-factors (~ 10 9 at 4K), which are important characteristics for electromechanical ex- 
periments at the quantum limit. We report the first cooldown of a bulk sapphire sample below 
superfluid liquid helium temperature (1.6K) to as low as 25mK. The electromagnetic properties 
were characterised at microwave frequencies, and we report the first observation of electromagneti- 
cally induced thermal bistability in whispering gallery modes due to the material T 3 dependence on 
thermal conductivity and the ultra-low dielectric loss tangent. We identify "magic temperatures" 
between 80 to 2100 mK , the lowest ever measured, at which the onset of bistability is suppressed and 
the frequency-temperature dependence is annulled. These phenomena at low temperatures make 
sapphire suitable for quantum metrology and ultra-stable clock applications, including the possible 
realization of the first quantum limited sapphire clock. 



Experiments to couple superconducting qubits based 
on Josephson junctions to microwave resonators (Q rs 
10 4 ) at cryogenic temperatures have been well repre- 
sented in recent scientific literature for a diverse range 
of circuit quantum electrodynamic applications. This in- 
cludes generating non-classical states of microwave cav- 
ities such as Fock states, where the limit on producing 
these non-classical fields is due to the finite photon life- 
time (or linewidth) of the resonator, as well as detecting 
a nanomcchanical resonator at or near the ground state 
[THZ] • Sapphire resonators are of particular interest for 
future experiments due to their extremely low loss, with 
electronic Q-factors of order 10 9 at 1.8K [SlU], and me- 
chanical Q-factors as high as 5 x 10 8 at 4.2K [TUl ITT] . 
The thermal, mechanical, and bulk electronic properties 
of sapphire have been characterised extensively over a 
wide range of temperatures from room temperature to su- 
perfluid liquid helium (1.6K) using WG mode techniques 
[121 113) . but have never been examined in the regime 
approaching the absolute zero of temperature. Sapphire 
resonators at millikelvin temperature have the potential 
to play an important role in the next generation of quan- 
tum electronics and metrology experiments by virtue of 
this anomalously high Q-factor. A significant body of 
research already exists in which sapphire has been used 
at cryogenic temperatures as a parametric transducer in 
an effort to reach the Standard Quantum Limit [TUl IT4T 
118] . Oscillators can be prepared in their quantum ground 
state due to very low thermal phonon occupation when 
T « hf/kg, where h and kfj are Planck's and Boltz- 
mann's constants respectively. For microwave oscillators 
such as those based on single-crystal sapphire resonators, 
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the corresponding temperature regime is in the experi- 
mentally accessible millikelvin range, making them ideal 
candidates for quantum measurement experiments. It 
is thus important to characterise such devices in this 
unexplored ultra-low temperature regime. In this Ar- 
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FIG. 1. Network analyser measurement of the WGIDo.o.o 
mode in transmission at 50 mK. The excitation power was 
varied in steps of 5dB from -50dBm to -15dBm, and the ob- 
served mode frequency was downshifted. The -50dBm and 
-45dBm curves are the highest in frequency and lie on top of 
one another. The threshold power is defined to be the power 
incident on the resonator which is sufficient to shift the mode 
frequency by one bandwidth from the 'unperturbed' lowest 
power measurement 

ticle we report on the first measurements of the elec- 
tromagnetic properties of a single-crystal sapphire res- 
onator at millikelvin temperature. The resonator used 
was a highest purity HEMEX-grade sapphire from Crys- 
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tal Systems, similar to that used in several Cryogenic 
Sapphire Oscillator (CSO) EH [20] and Whispering 
Gallery Maser Oscillator (WHIGMO) experiments [2H - 
[2^] . Here, we report on the observation of the low- 
est frequency-temperature turning points for Whisper- 
ing Gallery (WG) mode resonances ever measured, as 
well as making the first observation of a thermal bista- 
bility effect in sapphire for this ultra-low temperature 
regime. We give a model to predict thermal bistability 
threshold power and show that the effect is dependent 
on the thermal conductivity of the sapphire. Further- 
more, we show that the bistability effect may be sup- 
pressed by operating at a "magic temperature" , where 
a frequency-temperature turning point occurs. Thus, 
combining the low temperature operation with stabilisa- 
tion at a millikelvin frequency-temperature turning point 
gives the potential to realise a sapphire based frequency 
standard at the quantum limit (rather than the thermal 
limit which has been previously reported ;.22J). 

The sapphire resonator, a cylinder 5cm diameter x 
3cm height, was cleaned in acid and mounted in a silver 
plated copper cavity. The resonator is machined such 
that the anisotropy c-axis of the sapphire is aligned with 
the cylindrical z-axis. A radially oriented loop probe and 
axially oriented straight antenna were used to couple mi- 
crowave radiation in and out of the crystal. The cav- 
ity was attached to the mixing chamber of a dilution 
refrigerator with a copper mount and cooled to 25mK. 
The fundamental quasi-transverse magnetic Whispering 
Gallery modes WGH m ,o,o, with azimuthal mode num- 
ber m from 13 to 20, were characterised over a range of 
temperatures using a vector network analyser. We ob- 
served that particularly high-Q WG modes exhibited a 
hysteretic behaviour, which was thermal in nature. The 
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FIG. 2. The WGH2o,o,o mode in transmission for varying 
sweep time. Sweep direction is increasing in frequency in all 
cases. The linewidth of the mode narrows sharply with sweep 
time due to temperature dependence of permittivity of the 
sapphire. The governing equations for the lineshape are given 
in [25] 



frequency of the WG modes supported in the resonator 
are dependent on both the physical dimensions of the 
crystal and its permittivity, the latter effect being more 
than an order of magnitude stronger |11| . As the network 
analyser sweeps in frequency, heating occurs as power is 
deposited into the sapphire on resonance. The change 
in permittivity due to temperature causes a shift in the 
resonant frequency of the mode in the opposite direc- 
tion to the sweep. The result is an astoundingly nar- 
row, yet artificial linewidth with a sharp threshold. If 
the frequency was swept in the opposite direction, the 
mode frequency is shifted in the same direction as the 
sweep, and an artificially broadened linewidth would be 
observed. A similar effect in such dielectric resonators 
has only been observed at optical frequencies in fused 
silica microspheres [26 . It was shown that the "ther- 
mal bistability" caused either narrowing or broadening 
of the line resonance depending on the direction of the 
frequency sweep during measurement. Examples of opti- 
cal bistability are numerous in the literature [23 [23 [2E] , 
but are normally attributed to a x Kerr nonlinearity, 
which results in a threshold power for optical bistabil- 
ity that scales like Q~ 2 . Collot et al. [35] note that for 
mode Q-factors below 10 9 the thermal bistability effect 
dominates over the Kerr effect due to significantly lower 
threshold power. For quality factors in the range of 10 9 , 
the effects can be distinguished by the observed depen- 
dence of threshold power on Q. We find an excellent fit 
using a thermal model (see Eqn{3] and Figj4]) which has 
a threshold power that scales like <5 _1 , showing that the 
effect is clearly thermal in nature. A measurement of 
the WGH 2 o,o,o mode was made (see Fig. [T| which shows 
the first observation of this thermal bistability effect in 
a millikelvin sapphire resonator. A sharp threshold was 
observed, giving an FWHM linewidth of only 0.00173 Hz, 
which was strongly dependent on input power. Our ex- 
perimental apparatus was unable to sweep downwards in 
frequency, but the bistability effect is still observable by 
varying the sweep speed. Figure [2] shows the effect of the 
thermal bistability for a range of sweep times. Note that 
only the resonant peak moves; the longer sweep time re- 
sults in more time spent per measurement point, deposit- 
ing more power into the resonator and creating a larger 
apparent frequency shift and linewidth change. As suf- 
ficient heat is deposited into the resonator only on reso- 
nance, the off-resonance transmission does not depend on 
sweep time. It is possible to model the threshold power 
at which bistable behaviour becomes apparent. Consid- 
ering a temperature dependent fractional frequency shift 
of the WG mode of interest, Avjv = — aAT, where the 
temperature coefficient a is experimentally determined, 
we then expect thermal bistability for a threshold tem- 
perature rise of the resonator ATjj, = ~ ~ = ^=> K. An 
expression for the threshold power required to achieve 
this bistability is given by |26j : 

= C pP V eff AT th 
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FIG. 3. Temperature dependence of Q-factor (shaded circles) 
and frequency (empty circles) for the 20 th azimuthal order 
WG mode. Note that due to paramagnetic impurities in the 
sapphire, df/dT is annulled at 2.0897K. The Q-factor remains 
approximately constant over a wide range of temperature. 



where C p is the heat capacity of sapphire, p is its den- 
sity (4.0 g/cm 3 [H]), tt the characteristic heat diffusion 
time constant, and V e ff the effective volume occupied by 
the Whispering Gallery mode. The heat diffusion time 
constant in turn may be expressed as: 



T T = 



ImCp 
Ak 



(2) 



where I and A are the length and cross-sectional area 
of the sapphire segment, m is the mass of the sapphire, 
and k the thermal conductivity [30] . Finally, an expres- 
sion (independent of the heat capacity and thermal time 
constant) is derived by combining Equations [l] and [5] 

Ak m P ff , , 

p th = Tn- ?£L 3 
laLJ m 



Where m e ff is the mass of the effective volume occu- 
pied by the WG mode. The expression for the threshold 
power is now clearly only a function of the thermal co- 
efficient a, Q— factor, and thermal conductivity of the 
sapphire, as well as its dimensions. The thermal con- 
ductivity k was estimated by fitting an approximate T 3 
power law to extrapolate below IK from the data for sap- 
phire in Touloukian et al.[3T], giving k = 0.039T 28924 W 



TABLE I. Measured parameters for the thermal coefficient 
and Q 
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mains similar to that at liquid helium temperature be- 
cause the heat capacity follows a similar cubic law to the 
thermal conductivity, leaving the ratio unchanged (Eqn. 
[2J. However, the threshold for bistability is substantially 
lowered with respect to liquid helium temperature due to 
the reduction in thermal conductivity of the sapphire. 
To experimentally determine the thermal coefficient a, 
frequency measurements of several Whispering Gallery 
modes were made over a range of temperatures from 
25mK to 5.5K. The modes examined were fundamental 
quasi-transverse magnetic modes WGHi4, ,0j WGHig i0 ,o 
and WGH 2 o,o,o- The modes were excited using a vector 
network analyser at low power, typically —45 dBm. In 
this way, saturation of residual paramagnetic spins in the 
sapphire was avoided. The temperature of the resonator 
was controlled at a number of points between 25-5500mK 
using a Lakeshore Model 370 AC Resistance Bridge, and 
custom data acquisition software recorded the temper- 
ature, Q-factor, and frequency of the modes in trans- 
mission. The base temperature of the dilution refriger- 
ator was 23mK, and temperature control was stable to 
within several millikelvin. The temperature dependence 
of mode frequency was mapped to produce plots such as 
Fig. [3] and several temperatures were chosen to measure 
the threshold power at which thermal bistability became 
apparent. The temperatures were chosen to reflect a 
range of values for the thermal coefficient a, ranging from 
nominally zero near the frequency-temperature turnover 
point, to a maximum at the largest slope. As previously, 
we define the threshold power to be the power incident 
on the resonator which is sufficient to shift the mode 
frequency by one bandwidth from the unperturbed low 
power measurement. Equation[3]is then used to calculate 
the theoretical threshold power. The experimentally de- 
termined thermal coefficients are summarised in Table [TJ 
and a particular example of the measured and calculated 
threshold power is given in Fig. [4] for WGH2o.o,o- 

Clearly, operation at the frequency temperature turn- 
ing point is advantageous as the thermal coefficient a ap- 
proaches zero and the threshold power required for bista- 
bility approaches infinity. Temperature turning points 
have been observed with no bistability at input pow- 
ers up to 20 dBm from 4-9K in similar sapphire res- 
onators, and are caused by residual paramagnetic impu- 
rities such as Ti 3+ , Cr 3 +, Mo 3+ , V 3+ , Mn 3+ , and Ni 3+ 
present at concentrations of parts-per-billion to parts- 
per-million. The opposite sign effects of temperature- 
dependent Curie law paramagnetic susceptibility, and 
temperature dependence of permittivity [32 36 cause a 
turnover in the frequency-temperature dependence. Op- 
erating at this turning point or "magic temperature" al- 
lows frequency fluctuations due to temperature instabil- 
ity to be annulled to first order, and has been crucial to 
achieve state-of-the-art short term fractional frequency 
stability in CSOs in the past [9, 20, 32 . Our results are 
the first observation of temperature turning points below 
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FIG. 4. Predicted threshold power (using Eqn[3| as a func- 
tion of temperature for WGtUo.co- Shaded circles are the 
measured threshold values for 200, 630, and 2100 mK. 



the boiling point of liquid helium. Table[II| lists the magic 
temperatures (turning points) measured for a range of 
WG modes. Operation of the WHIGMO at a millikelvin 
magic temperature rather than its current «8K would 
have the benefit of reduced thermal noise floor of the 
maser, with potential operation at the quantum limit. 
Operation at a millikelvin turning point, where the ther- 
mal coefficient a is vanishing, would minimise the effects 
from the considerable heating due to the large (>10 dBm) 
input power required to saturate the pump transition of 
the maser. This is similarly advantageous for CSOs [2U] 
which typically circulate large amounts of power through 
the sapphire resonator, at least several milliwatt, while 
the cooling power at base temperature of a dilution re- 
frigerator is only on the order of several hundred fxW. 
Additionally, high power operation of quantum limited 
transducers could be attained at these temperatures. 



TABLE II. Measured 'magic temperatures' for a range of 
quasi-transverse magnetic WG modes. Several modes close 
to the Fe 3+ centre frequency exhibited strong distortion and 
could not be accurately tracked to determine the turning 
point. The m = 17 mode required large power to excite and 
could not be measured below 80mK due to heating effects. 



than an order of magnitude lower than previously 
achieved, the first measurement of thermal bistability in a 
microwave sapphire resonator at these temperatures, and 
the first observation of millikelvin frequency-temperature 
turning points. We give a model for the thermal bista- 
bility threshold power, and show that it is closely de- 
pendent on the material properties of the sapphire res- 
onator. We propose several reasons the effect has not 
been previously observed in this system. Typically the 
CSO/ WHIGMO is operated very close to a frequency- 
temperature turning point where the temperature coef- 
ficient a is small, and the threshold power for bistabil- 
ity becomes significantly larger than normal operational 
power levels. In the experiments reported in this paper, 
the resonator was operated at low temperature, far from 
a turning point where the temperature coefficient was 
large, leading to a lower and readily observable thresh- 
old power. Frequency-temperature turning points as low 
as tens of millikelvin were observed for WG modes in 
the resonator, and we show that the thermal bistability 
effect can be suppressed by operating at these "magic 
temperatures" . Additionally, mode Q-factors remained 
high and were comparable to their usual values at higher 
temperature, ruling out the existence of extra loss mech- 
anisms in sapphire in the millikelvin regime. We con- 
clude that single-crystal sapphire is an excellent candi- 
date for both clock applications, and quantum metrology 
of macroscopic systems at temperatures approaching ab- 
solute zero. 
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In summary, we report the first characterisation of a 
single-crystal sapphire resonator at temperatures more 
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